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We carried out a finite-size scaling analysis of the restricted solid-on-solid version of a recently 
introduced growth model that exhibits a roughening transition accompanied by spontaneous symme- 
try breaking. The dynamic critical exponent of the model was calculated and found to be consistent 
with the universality class of the directed percolation process in a symmetry-broken phase with a 
crossover to Kardar-Parisi-Zhang behavior in a rough phase. The order parameter of the roughen- 
ing transition together with the string order parameter was calculated, and we found that the flat, 
gapped phase is disordered with an antiferromagnetic spin-fluid structure of kinks, although strongly 
dominated by the completely flat configuration without kinks. A possible interesting extension of 
the model is mentioned. 
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I. INTRODUCTION 

It has been realized that nonequilibrium interacting particle systems are capable of exhibiting very interesting 
and unusual phenomena in (1+1) dimensions, such as single-defect and boundary induced phase transitions [Q-^. 
Moreover, it is known that some of these phase transitions are accompanied by spontaneous symmetry breaking 
(SSB) in which some macroscopic observable of the model behaves in the steady state asymmetrically with 

respect to what it would be expected from the microscopic rules governing its dynamics. All this is rather unusual 
for equilibrium one-dimensional systems of particles interacting through short range forces only, and there have been 
serious attempts towards the understanding of these phenomena, in particular that of SSB jsj. It appears that among 
the many ingredients favoring SSB to take place, unbounded noise is a key one, whether it comes from the introduction 
of defects, boundary terms or any microscopic rules rendering a system lacking in detailed balance. 

Growth models provide a suitable theoretical framework for the investigation of a gamut of model systems ||^ , and 
it turns out that they can be used to investigate the above mentioned phenomena as well. In a recent work 0, a 
class of growth models addressing both the questions of a roughening transition in one dimension and of SSB was 
introduced. Concerning SSB, these authors have inquired about the necessary ingredients in a model in order for 
it to show SSB, and they found it possible to have SSB characterized by a nonconserved order parameter in a ring 
geometry, as opposed to the situation in a related model showing the same characteristics but now with the SSB 
associated with a conserved (i.e., slow) order parameter and in the presence of conspicuous boundary terms 
In the unrestricted version, the models can be related to a directed percolation process thus sharing its exponents, 
a fact that was confirmed by Monte Carlo simulations 0. In the restricted solid-on-solid (RSOS) version no such 
relationship exists. It is possible, however, to arrive through a site-link transformation at an equivalent driven diffusive 
system in which two types of oppositely charged particles diffuse asymmetrically and are continuously created and 
annihilated in pairs. 

In this work we have proceeded to a further investigation of the RSOS version of the models first proposed in [Q. 
Our study is based on the mapping of the master equation governing the dynamics of the associated reaction-diffusion 
process into an imaginary-time Schrodinger equation, the Hamiltonian of which is that of a spin S = \ non-Hermitian 
quantum chain ||^. This allows us to employ standard finite-size scaling (FSS) techniques to the resulting stochastic 
process, in the same way as one is used to do with Hamiltonian theories or with the transfer matrices of classical 
^ . spins systems |p. In this way we were able to analyze the time evolution operator for chains of sizes up to 16 sites, 
- - ' calculating their spectra and stationary states. The paper is organized as follows. In Sec. II we present the basic 
formalism in which we work, comment on its physical contents and derive the time evolution operator for the particular 
process we are interested in. In Sec. Ill we show and discuss the finite-size data we obtained for the dynamic critical 
exponent, the order parameter of the roughening transition, and the string order parameter. Finally, in Sec. IV we 
conclude and indicate some directions for further investigation. 
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II. THE TIME EVOLUTION OPERATOR 



We are going to focus on the stochastic particle system associated with the growth model, more specifically on 
its stochastic transition or intensity matrix. There is quite a number of different ways of writing down the master 
equation in operator form, some more suited to the study of symmetries some others envisaging a perturbative 
approach . Here we give a brief derivation that parallels that of |^ . 

Let us attach to each site ^ of a one-dimensional lattice A C Z of volume |A| = L a stochastic variable ne taking 
values in the set of states w = {0, 1, . . . , iV — 1}. Denoting by P(n, t) the probability of realization of a particular 
configuration n = (ni, ?i2, . . . , n^) G O = uj^ at instant t, we write the master equation as 



dPjn, t) 
dt 



nen 



[r(n,fi)P(n,i)-r(fi,n)P(n,i)] 



(1) 



where r(n, n) > is the rate for the transition n ^ fi. When only binary collisions intervene, we write r(n, n) 
for the elementary process (a, b) (c, d), and the master equation reads 
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where the primes indicate that the free channel (a, b) = (0, 0) = (c, d) should not be considered in the summations and 
the additions in the indices of the rates and in the arguments of P(n, t) are all taken modulo N. Periodic boundary 
conditions on A will be understood in what follows. 

We now introduce vector spaces in the description of Eq. To do this we turn uj = {0, 1, . . . , iV — 1} into uj — 
and n into | n) = | ni , 712 , . . . , ni, } ^ ft 
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. Taking the orthonormal basis {|n)} for f2 we write 
|P(i)) = 5]P(n,i)|n) 



(3) 



for the generating vector of the probabilities P{n, t) — {n\P{t)). We are in this way providing the space of generating 
functions with a Hilbert space structure. Next we endow the space of linear transformations of uj with the canonical 
basis of matrices with elements {E°-^)ij = SaiSbj, < a, 6, z, j < iV — 1, and introduce the operator X such that 
Xij = 6i+ij, = 1. A little reflection and comparison show that within this settings we can write Eq. (||) in the 
form 
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H is but the infinitesimal generator of the Markov semigroup U{t) = cxp{—tH) of the continuous-time Markov chain 
defined by the set of rates F^^- From Eq. (|^) we see that H is a N'" x , usually nonsymmetric, and very sparse 
matrix. One says that Eq. (0) is an imaginary-time Schrodinger equation, the fact that H has not the physical 
significance of energy and in general observes iJ^ 7^ H notwithstanding. It is however useful and intuitive to think of 
H as a Hamiltonian whose quantum fluctuations govern the time fluctuations of the classical system of particles. 

The process we are interested in is a growth model defined as follows . Let hi G N he the height of a surface at 
site ^ £ A. The surface evolves by attempting, sequentially and at randomly chosen sites, adsorption of an adatom 
he ^ hi + 1 with probability qdt, and desorption of an adatom he — > imn{he-i, he} or he —^ imn{he, he+i} each with 
probability ^(1 — q)dt. We now impose the RSOS condition \he+i — he\ < 1,V^ G A, which suggests the use of link 
variables ce = — he E {—1,0, 1}. In order not to overload the notation with unnecessary pluses and minuses, 
let us take the values of ce modulo 3, mapping q = —1,0, 1 into Ci = 2,0, 1, respectively. In this representation the 
growth process can be described by the set of transition rates 
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and = g. 



According to Eq. (J5|) we write for this process the time evolution operator as 
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with the two-body stochastic transition matrix given exphcitly by 
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where we have ordered the two-site basis vectors antilexicographicaUy, i.e., |0,0) |1, 0) ^ . . . -< |1,2) -< |2,2), and 
the dots indicate nuU entries. 

Before proceeding to the next section, it is worth mentioning some properties of H. As a stochastic transition 
matrix, conservation of probabihty flux require s Hij = 0, which in turn constraint the diagonal elements to be 
given by Hjj = — J2i=£j ^ij j compare with Eq. (|^ The first of these conditions imply the existence of a trivial left 
eigenstate with zero eigenvalue, 



(9) 



with (rjjri) = the cardinality of the state space. This special vector plays a role in the determination of expectation 
values, for one can write the average of an observable A{n) with respect to the probabilities P(n, t) with the aid of 
|f2) as 



(A)(i) = EA(n)P(n,t) = (l]|A|P(t)). 



(10) 
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We expect physical observables of the classical system of particles to be diagonal in the natural basis {|n)}, once 
they all have to commute. H is obviously not diagonal in this basis, it is not an observable of the system. Eq. ( [l0| ) 
summarizes an important difference between quantum physics and the kind of classical physics we have here: it is 
that expectation values are linear in \P{t)), not bilinear, the elements of \P{t)) being probabilities themselves, not 
probability amplitudes. 

Besides these general properties, our H in Eq. (|^) is additionally translation invariant, due to the homogeneous 
rates in a ring geometry, and possesses a U(l) symmetry labeled by the total charge Q = Q^+^ — Q^^-*, which is 
conserved along the process. These symmetries allow us to block-diagonalize H and to write it as the direct sum 
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= -L fc=0 



where Q and k label the U(l) and momentum eigensectors, respectively. Each U(l) sector represents a closed class 
of the stochastic process, the corresponding block matrix H'^ being itself a stochastic transition matrix governing 
the dynamics within the given sector. We say that our process is decomposable, nonergodic, and that the U(l) label 
classifies its 2L + 1 irreducible, closed classes. The momentum label fc, however, is introduced here solely in order 
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to take advantage, numerically, of the further reduction of order 1/L on the sizes of the blocks it furnishes, the 
physically relevant momentum sector being the one with k = 0, since the zero eigenvalue together with the completely 
flat surface of the model is in this sector. One can look for a relation E{k) oc for the real part of the eigenspectrum, 
but we do not expect to extract useful information, e.g., about 9, from this relation for small lattice sizes. Moreover, 
the interpretation of such a relation for the low lying excitations of the process as a dispersion relation for quasiparticles 
of a free field theory would be somewhat cavalier; see 



iiT|. 



III. FINITE-SIZE SCALING 

According to Eqs. (^), as q increases creation of H — as well as annihilation of — h pairs increases while the 
remaining processes induce segregation of particles. This corresponds to an increase in adsorption and in the growth 
rate of islands, leading to rougher configurations. As q lowers, increased annihilation of H — pairs together with more 
symmetric diffusion of O's flatten the surface; see Sec. Ill C. Around the critical point q — qc separating these two 
phases the correlation lengths of the infinite system behave like 

(X el « \q-qcr« cc \q~qcr^'. (12) 

Here is the correlation length in the time direction while e_L is the one in the spatial direction, with and i'± the 
corresponding critical exponents and = the dynamic critical exponent. 

For finite systems of size L, according to the usual FSS assumptions we expect the scaling 

Cllx^i'S (13) 

to hold when q — qc,L, the finite version of the critical point gc, with 9l the finite version of 9. On general grounds 
one expects limL^oo qc.L = qc and lim^^oo 9^ = 9. Therefore from Eqs. ( p^ ) and ( p^ ) we obtain the relations 

In [e||,L(gc,L)/e[|,L-l(gc,L)] ^ In [ClL + liqc,L)/^lLi<lc,L)] ^ , . 

ln(L/L-l) ^ ln{L + l/L) ^ ^' ^ ^ 

which through the comparison of three successive system sizes furnishes simultaneously gc.L and 9^. Of course 
A-i = Re{A£;L}, with AEl = ^ the first gap of since for stochastic transition matrices E^^^ = by construction. 



A. Dynamic critical exponent 

We have calculated the gaps AEl in the Q = 0, fc = sector through the use of the Arnoldi algorithm |]l5| , |l6| . 
This is a Krylov subspace projection technique that effects the reduction of a general nonsymmetric matrix to upper 
Hessenberg form, the eigenpairs of which converge variationally to that of the original matrix with the algebraically 
larger part of the spectrum converging first. In order to save memory we have used a restarted version of the algorithm, 
in which we fix the dimension of the Krylov subspace and use some of the approximate eigenvectors obtained in one 
iteration as the starting vectors for the next iteration, until convergence is obtained to the desired accuracy. In this 
way we were able to handle matrices of orders up to 324 862 with up to ~8 x 10^ nonzero entries keeping the Krylov 
subspace always with fewer than 64 vectors. 

The results we have obtained are summarized in Tabs. || and |l|. The extrapolated values in the last line of these 
tables were obtained through the Bulirsch-Stoer extrapolation scheme with lubst the free parameter of the 
algorithm chosen over a certain range so as to optimize the converge of the finite-size data. 

In applying Eqs. ( |l^ we found two consistent, converging sets of data, the first one, shown in Table | and marked 
with a prime, realizing the first equality in Eqs. (|l^) only approximately, and the second one, that in Table || and 
marked with two primes, realizing it exactly. This behavior is illustrated in Fig. |^. 

The first set of data exhibited a rather smooth convergence in both the values of qc.L and 9l, while the second set 
behaved more irregularly. The values in Table | indicate a second order transition taking place around q'^ ~ 0.1875 
with a dynamic exponent of 9' ~ 1.585. This value of 9' is compatible with the exponent of the directed percolation 
process, for which the most accurate value to date, obtained by Monte Carlo simulations, seems to be 9^? — 1.580 75(3) 
Ipsf; it should be mentioned that in an evaluation of 0dp more closely related to ours it has been found a value of 
I, and that it is not clear why such discrepancies arise in the value of 0dp calculated by different 
methods. 
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The interpretation of the data in Table y is touchier. We can see that the values of g"^ converge at a reasonable 
rate to the extrapolated limit ~ 0.1932, which is different from the previously found q'^ ~ 0.1875. We believe 
however that if we have had access to larger lattice sizes, we would have observed — > g^, since the values of q'^ 
vary less. We thus trust the value of q'^ = 0.1875(1) as our best estimate for the critical point. The situation with 
the critical exponent 9'^ is different: it seems to be converging to a completely different value than 0dp- In fact, this 
behavior was to be expected, for it has been found that in the rough phase q > qc the exponents of the process 
are those of the Kardar-Parisi-Zhang universality class ||2^, in particular Okpz = f ■ As can be seen from Table ||, 
the first few values of 9'1 show a monotonic increasing behavior up to 9g ~ 1.4649, but then the sequence begins to 
decrease to reach the bottom value of ~ 1.4426. A partial extrapolation of the first four points, 6 < L < 9, gives 
9" = 1.47, while a partial extrapolation of the rest of the points, 10 < L < 15, gives 9" = 1.42. This absence of 
monotonicity of the finite-size data for 9'1 is quite unusual, and we have not yet a clear clue to this behavior. Even 
so, we can take the set of values of 9'1 in Table |l| as indicating the presence of a critical region for q > qc with an 
exponent 9" ^ 9bp, possibly with 9" = 9kpz- 

Concerning the exponent I'x, we found it not possible to apply the standard approach p^Jl^ to obtain its value 
because the derivatives of RejAi?!,} with respect to q evaluated at the points qc,L change sign for some pairs of 
lattice sizes, thus preventing us from taking logarithms. The finite-size sequences obtained with the absolute values 
of these derivatives as well as with those obtained with the derivatives of the absolute values of the gaps also failed 
to converge, so that we were not able to obtain an estimate of from our diagonalizations. 



B. Spontaneous symmetry breaking 

Given that the surface suffers a transition from a flat phase to a rough phase, it is natural to think of an order 
parameter which measures this transition. Moreover, it is interesting to have an order parameter taking into account 
the symmetries of the process, which are besides the translational and U(l) symmetries, a Zoo symmetry related to 
the fact that the microscopic dynamics is invariant under an arbitrary integer shift hi he + n in the heights. A 
proper order parameter is given by 0| 

Mi = lX](-l)''' (15) 

£=1 

which is fast in the sense of not being conserved by the dynamics. The choice of such an order parameter anticipates 
the interpretation of the roughening transition (actually, of the flattening transition) as the result of a spontaneous 
break of the Zoo symmetry, for while in the rough phase all heights are exploited evenly, in the flat phase the system 
spontaneously selects one fiducial level around which the heights fluctuate. One then expects Ml to be finite in the 
flat phase while vanishing in the rough phase due to canceling fluctuations. 

We have calculated Ml for even L between L = 8 and L = 16; our finite-size data together with the extrapolated 
values appear in Fig. ^. From that figure we clearly see the transition taking place around q — 0.190, although the 
precise determination of the critical point is not possible from this figure. We have not found the signature of two 
different qcS in our data for M — lim^^oo Ml, which we regard as an indication that should indeed tend to q'^ as 
L ^ oo. 

The order parameter M vanishes around 9 ;$ <Zc as M oc (qc — q)^. The plot of InM x \n{qc — q) for the points 
0.12 < q < 0.18 of the extrapolated curve in Fig. ^ together with a linear regression fit appear in Fig. |. We found 
an exponent 77 = 0.57 ± 0.03, which compares well with previous results in the literature: in the first of the papers by 
Alon et al. 0, 77 has been evaluated as 0.55 ± 0.05, while in a recent simulation of a model of yeastlike growth of fungi 
colonies with parallel dynamics it has been found that ~ 0.50 [^. Also in a certain line in the phase diagram of 
a one-dimensional next-nearest-neighbor asymmetric exclusion process closely related to these growth models it has 
been found that 77 = 0.54 ±0.04 In the second of the papers by Alon et al. 0, however, the more accurate value 
rj = 0.66 ± 0.06 has been published, pushing the estimate to a somewhat higher value. Recent preliminary Monte 
Carlo simulations of ours, on the other hand, suggest a typical Ginzburg-Landau scenario for the symmetry break in 
these models, which would then predict an 77 = i. We believe that more extensive simulations can settle this point, 
and work is being done in this direction. 



C. String order parameter 

In order to better understand the nature of the roughening transition, let us look at some typical microscopic 
configurations of the model. The roughest possible surface is given in the links representation by | + + • • • H ), 
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the state of maximal height in the Q = sector of the dynamics. In spin language, this state corresponds to two 
domains separated by two antiferromagnetic (AF) kinks. The second roughest possible surface configurations are 
given by those with a pair of O's, e.g., | + ■ • • + — ■ • • — ) or |0 + ■ • • H — ■ • • — 0). From this example and the rates 
in Eqs. ^ it becomes clear that the flattening process is induced by AF kink annihilation, while diffusion of O's 
introduce surface shape fluctuations. It is important to notice that absence of desorption from the middle of smooth 
terraces enforces a certain AF order among the particles, for pairs are created only as H — pairs, never as — h pairs, 
and since there are no H — ^ — h reactions (which would violate the RSOS condition), we see that this order persists 
as long as pairs survive annihilation. In the rough, high q phase the + particle will preferably move leftwards, while 
the — particle will prefer to move to the right, eventually leaving a pair of O's in between which then generates another 
H — pair, thus leading to rough configurations like, e.g., | • • ■ + +0 • • • — — • • ■). In the fiat, low q phase difi^usion of 
particles becomes more symmetric and segregation of particles less efficient, and we thus expect to observe a more 
uniform distribution of H — pairs along the lattice than in the high q phase. The completely fiat surface without any 
AF kink is only attained at g = 0. 

In the context of two-dimensional RSOS crystal growth models and the Haldane conjecture |^^, it has been 
predicted ||2^ and subsequently extensively verified |2^-[30|] that a particular type of long-range order exists in the 
spin S = 1 antiferromagnetic isotropic Heisenberg (AFH) chain. In the ground state of the AFH chain, this order 
may be viewed as made up of (not necessarily closely) bound H — dipoles interspersed among the O's, forming what 
has become known as an AF spin fluid. The order parameter that identifies this type of order is the so-called string 
order parameter |24l , 



0^(^)= 5^exp 



ITT y 



(16) 



where the brackets indicate expectation value in the ground state. In the gapped, Haldane phase of the AFH chain 
as well as in the disordered flat phase of the models in one has lim^^oo 0%[i) ^ 0. 

From what we said above it is clear that the string order is just the kind of order we expect to observe in the 
fiat phase of our model. We have thus calculated the steady state expectation value of the string order parameter 
( |l6|) with S^-^ = c„ and i = L/2 + 1, the maximum distance in a ring geometry, and the results appear in Fig. IJ. 
Our extrapolations did not perform well for this set of data, and are not showed in this figure. The general trend, 
however, is quite clear: above the critical point, q > qc, 0^{£) strongly tends to zero, while for q < qc we have 
lim^^oo 0^{£) ^ 0. The point g = is special, for at g = the completely fiat surface becomes an absorbing state and 
0^{£) = exactly at this point. As q grows from zero, the AF spin fiuid begins to form and 0^{£) grows accordingly, 
until a,t q — qc the asymmetry in the diffusion rates for the particles disrupts this AF spin fluid structure, ordered 
domains begin to prevail and the string order vanishes. We then see that the string order parameter clearly reveals 
the mechanism of the roughening transition as the unbounding of the fluid antiferromagnetic pairs in favor of the 
formation of ordered domains. 

IV. SUMMARY AND CONCLUSIONS 

In summary, we carried out a finite-size scaling study of the roughening transition in a class of one-dimensional 
RSOS models which also presents spontaneous symmetry breaking. We found that at the critical point qc — 0.1875(1) 
the transition occurs with a dynamic exponent compatible with that of the directed site percolation process, for 
which we have the estimate 6 — 1.585(1), and that above qc there is a critical rough phase most probably with KPZ 
exponents. Unfortunately, we were unable to calculate a second exponent from our diagonalizations. This might be 
due to the nonhermiticity of the operator H, which might have caused an unusual nonmonotonicity in the values of 
the gaps with the parameter q, thus preventing us from obtaining i'±. This lack of monotonicity has already been 
reported in the literature where it has also been noticed the slow convergence of the finite-size data towards 

the infinite volume limit. 

The order parameter M was found to vanish like M oc (qc — q)'' for q qc with an exponent 7] = 0.57 ± 0.03, in 
agreement with previously found values of rj [0,^ 22 1. The calculation of the string order parameter revealed that 



the fiat, gapped phase of the model is a disordered phase analogous to a Haldane phase, with the stationary state 
presenting an antiferromagnetic spin-fiuid structure of kinks, although dominated by the completely fiat surface with 
no such a structure. The roughening transition may thus be understood in the links representation as the unbounding 
of the fiuid antiferromagnetic pairs in favor of the formation of ordered domains, which then begin to blend together 
providing the surface with a finite growth velocity. 

It is possible to push further the investigation of this class of models in one definite way. The idea is to allow for 
an explicit break of the symmetry in the set of rates Eq. (ph by the following artifact [p2|. In the particle scenario. 
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we double the number of sites, introducing between two successive links q, q+i a noninteracting flag variable m^^i 

taking two possible values, call them + and — . This variable will mimic the pseudospin (—1)''*. We then allow the 
rates of our modified model to be parametrized by, besides q, a chiral symmetry breaking field u € [—1,1] such that 
now the rates depend on the quantities 

p = (1 — u)q and p — (I + u)q 

according to F"^^ — r°^(p) and r°^^ ~ ^triip)^ roles of p and p interchanged whenever rj^° ^ for a given 
Tf^ ^ 0. For example, rjt° = \{l - p) and rj^° = ^(1 - p), but T^tl ^ p and T°^l = p. The choice of which 
combination of values of to's in the new rates will pick a p or a p with respect to the original rates is immaterial, for 
letting u —^ —u exchange their roles. For this process one may look at the order parameter 

1 ^ 

£=1 

to see whether one finds a spontaneously symmetry-broken phase. It may happens that for some values of the field u 
one gets spinodal points, and that these points are associated with unusual dynamic exponents, e.g. 0=1, once they 
have already appeared in one dimensional driven diffusive systems [0 . This problem is currently under investigation. 
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TABLE L Finite-size data for the region where a minimum is observed. The numbers between parentheses represent the 
estimated errors in the last digit of the data. 



L-1,L,L + 1 


Q'c.L 


O'l 


5,6,7 


0.190 461(1) 


1.794411(1) 


6,7,8 


0.190 676(1) 


1.749 008(1) 


7,8,9 


0.189 988(1) 


1.718 350(1) 


8,9,10 


0.189 294(1) 


1.695 320(1) 


9,10,11 


0.188 753(1) 


1.677170(1) 


10,11,12 


0.188 355(1) 


1.662 589(1) 


11,12,13 


0.188 076(1) 


1.650 766(1) 


12,13,14 


0.187 885(1) 


1.641 112(1) 


13,14,15 


0.187761(1) 


1.633 190(2) 


14,15,16 


0.187 683(1) 


1.626 655(2) 


Extrapolated 


0.1875(1) 


1.585(1) 


[wbst] 


[7.329] 


[2.543] 


TABLE II. Finite-size data for the region where crossing occurs. The data without an associated error are believed to be 


correct to the figures shown. 




L-1,L,L + 1 


// 

<lc,L 


d'l 


5,6,7 


0.295 749 


1.445 967 


6,7,8 


0.275 660 


1.458 235 


7,8,9 


0.261781 


1.463 535 


8,9,10 


0.251 648 


1.464 937 


9,10,11 


0.244 053 


1.463 830 


10,11,12 


0.238 186 


1.461059 


11,12,13 


0.233 511 


1.457203 


12,13,14 


0.229 689 


1.452 672 


13,14,15 


0.226 497 


1.447 753 


14,15,16 


0.223 784 


1.442 647 


Extrapolated 


0.1932(1) 


see text 


[wbst] 


[1.815] 


see text 
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FIG. 1. Variation of the exponent 9 with q. The values of 0^^^ and 9^'^^ were obtained from the first and second expressions 
in Eqs. respectively, using the triplet of lengths L — 1, L, L + 1 = 13, 14, 15. The inset shows that the difference between 

their values reaches a minimum around q = 0.188 and vanishes around q = 0.226. The finite-size sequences obtained from both 
the minima and the crossings converge to well defined limit values, see Tables | and |]lj and the text. 




FIG. 2. Order parameter Ml for even L between L — 8 and L = 16 together with the extrapolated curve. 
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